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We present an experimental and theoretical study of the phase behavior of a binary mixture of 
colloids with opposite adsorption preferences in a critical solvent. As a result of the attractive and 
repulsive critical Casimir forces, the critical fluctuations of the solvent lead to a further critical point 
in the colloidal system, i.e. to a critical colloidal-liquid-colloidal-liquid demixing phase transition 
which is controlled by the solvent temperature. Our experimental findings are in good agreement 
with calculations based on a simple approximation for the free energy of the system. 

PACS numbers: 82.70.Dd, 68.35.Rh, 81.16.Dn 



The structure, phase behavior and dynamics of col- 
loidal suspensions is largely determined by the interac- 
tions between the colloids and considerable effort goes 
into controlling the pair potentials in such systems. This 
can be achieved, e.g. by adding salt to the solution in 
order to screen electrostatic interactions pQ, by adsorb- 
ing charged or polymeric species onto the surfaces of the 
colloids PP, or by adding a depletion agent to the so- 
lution, such as a non-adsorbing polymer [TJ [2]. While 
these methods have been successfully employed in many 
experiments, they can not be used when colloidal inter- 
actions have to be changed in-situ and in a reversible 
manner. In such situations, external fields may be used 
and have the additional advantage of allowing for time- 
dependent control of the colloidal interactions. Exam- 
ples are superparamagnetic or dielectric colloidal parti- 
cles, where the structure and phase transitions can be 
easily controlled by external magnetic and electric fields 
[31 13] or temperature-responsive microgel colloidal par- 
ticles, which change their size when the temperature is 
varied by several degrees [5] . 

Only recently, critical Casimir forces have been demon- 
strated to provide an alternative approach for the in-situ 
control of colloidal forces [5HH]- Such interactions arise 
between particles that are immersed in a fluid mixture 
which is near to its critical point |10) . The confinement 
of the critical concentration fluctuations between the sur- 
faces of a pair of colloids leads to a force whose range and 
amplitude strongly depend on the correlation length in 
the mixture and thus shows a strong temperature de- 
pendence [TTJ Q2]. In addition, critical Casimir forces 
can be changed from attractive to repulsive by changing 
the adsorption preference of the colloidal particles for the 
mixture's components, which determines the correspond- 
ing boundary conditions (BC) 8 . So far, all experiments 
have been performed with colloidal systems where all par- 
ticles have the same BC which then leads to particle ag- 
gregation or glass formation in concentrated suspensions 

[BUM!]. 

In this Letter, we investigate the phase behavior of a bi- 
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FIG. 1: Critical Casimir potentials obtained from Eq. 
with the scaling functions from Ref. [6], for two 3 /xm sized 

particles with symmetric (H — h, ) and antisymmetric (H — ) 

boundary conditions at various different temperatures. Inset: 
a) Experimental snapshot of the demixing process of the crit- 
ical water - 2,6-lutidine mixture, b) Snapshot of the binary 
colloidal mixture, with one species labeled with a fluorescent 
dye. The scale bar corresponds to 3 (im. 



nary colloidal mixture with opposite boundary conditions 
that are dispersed in a binary critical mixture. Due to the 
attractive and repulsive critical Casimir forces for parti- 
cles with the same and opposite BC respectively, we ob- 
serve a critical colloidal-liquid-colloidal-liquid (CL-CL) 
demixing phase transition which is entirely driven by the 
critical concentration fluctuations of the solvent; i.e. the 
critical point in the solvent mixture generates the phase 
transition in the colloidal system. This demonstrates a 
strong coupling of different critical phenomena over more 
than two decades in length scales. Our measurements are 
in agreement with the theoretically calculated phase dia- 
gram obtained from a simple approximation for the free 
energy of the system. 

The critical Casimir potential between two colloidal 
spheres with radius R with center-center distance r is 
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FIG. 2: Snapshots of the particles in binary colloidal mixtures of two different compositions: in panels a-e, xa = 0.54 and in 
f-j, xa = 0.29. The temperature is approaching T c , as one goes from a to e and also from f to j, respectively. The last image 
in each column is a snapshot of the system at a longer time. Panels a-d are obtained after lh at the temperatures AT = IK, 
0.40K, 0.08K and 0.01K and e is obtained after 20h at AT =0.01K. Panels f-i are obtained after lh at AT = IK, 0.09K, 0.05K 
and 0.02K and j is obtained after 22h at AT =0.02K. 



given by [BJ: 



where z = r — 2R is the distance between the particle's 
surfaces and £ = £ (AT/T C )~ - 63 is the solvent bulk cor- 
relation length of the binary mixture with a critical tem- 
perature T c . The value £o is a microscopic length scale 
which characterizes the range of the molecular interac- 
tions in the mixture. The function #.y(z/£) is a universal 
scaling function with (i = j) corresponding to symmetric 
and (i ^ j) antisymmetric BC. Note the factor 1/2 in Eq. 
([IJ, since we used the scaling function for a sphere- wall 
geometry [12"] . 

Our experimental system is composed of a binary mix- 
ture of silica spheres with slightly different diameters 
a a = 2.8 /jm and ob = 3.4 /im which are suspended in 
a water-2,6-lutidine mixture at critical composition (lu- 
tidine mass fraction of c c L = 0.286). The solvent is char- 
acterized by a lower critical temperature T c = 307K (i.e. 
the mixture is in its homogeneous state below T c ) and 
£o ~ 0.2nm±0.02nm, as determined from light scattering 
experiments [16]. Initially, both species of particles (de- 
noted A and B) have a strong adsorption preference for 
water, i.e. with (— ) BC, due to their high surface charge 
[17] . In order to reverse the BC of the A- type particles, 
they are functionalized with silane according to a proce- 
dure described in |21j . This treatment renders the parti- 
cles hydrophobic, i.e. (+) BC, as inferred from the fact 
that they preferentially stay in the lutidine-rich phase of 
the demixed binary solvent. The entire colloidal suspen- 
sion is contained in a temperature-controlled thin sample 
cell in which all the colloids sediment towards the bottom 
plate and form a rather dense monolayer. The sample 
plane is imaged onto a camera and the particle positions 
are determined with video microscopy. To facilitate the 
discrimination of particles with opposite BC, the B-type 



particles are labeled with a fluorescent dye (see inset of 
Fig.JT]). We should mention that for the particle size ra- 
tio used here, no demixing due to depletion interactions 
has been observed. We should also mention that critical 
Casimir forces also arise between the particles and the 
bottom of the container. However, these do not influence 
the observed demixing because these forces only change 
the average particle height above the substrate and so, 
e.g., change the particle mobilities due to hydrodynamic 
interactions [HI [19] . In case of the silica particles used 
here, such variations in hight can not be experimentally 
resolved due to the large gravitational forces which lo- 
calize the particles at roughly the same hight above the 
surface. As a result, the change in the mobility of the 
particles is negligible. Note also that bridging [6] [20] 
between the particles and the substrate does not occur 
because we use a mixture with a critical composition and 
always stay below the critical temperature. 

Fig. [I] shows the critical Casimir interaction <&ij(r) for 
different temperatures between particles with 3 fim di- 
ameter and with symmetric (H — h, ) and antisymmet- 
ric (H — ) BC as calculated from Eq. 0. Additionally, 
electrostatic contributions have to be considered, which 
depend on the particle's surface charge and the Debye 
screening length of the solvent with k~ x w 12nm for 
a critical water-lutidine mixture [BJ. Because the posi- 
tion of the first peak in the pair correlation functions of 
the colloidal mixture are independent of the particle den- 
sity in the temperature range covered in our experiments, 
the electrostatic interactions are well approximated by a 
hard sphere potential with an effective particle diameter 
a e ff = <7i + 0.1/Ltm (dashed vertical line in Fig. [I]). 

Fig. [2] shows temperature dependent configurations of 
a binary colloidal system with total density pa&a 2 + 
Pb&b 2 = 0.70 ± 0.05, where pa and pb are the (two- 
dimensional) surface number densities of the A (or- 
ange) and the B (black) particles for two mixing ratios 
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xa = Pa/{pa + Pb) = 0.54 (Figs. [2]a-e) and x A = 0.28 
(Figs. [2]f-j). Each temperature change is followed by a 
waiting time of about lh to allow the particles to rear- 
range accordingly. When increasing T from room tem- 
perature towards T c , up to AT = OAK, the particles re- 
main homogeneously distributed for both compositions 
(Figs.[2]a,b). This clearly suggests that in this tempera- 
ture range, the critical Casimir pair interaction energy is 
less than or equal to the thermal energy as also confirmed 
by Fig. [I] Further approaching the critical point, the crit- 
ical Casimir forces increase and we observe large struc- 
tural changes in the colloidal mixture. For xa — 0.54 the 
system changes from a random distribution to an almost 
bicontinuous network. We see this beginning to set in, in 
Fig.[2]c, which corresponds to AT = 0.08K . Further in- 
creasing the temperature of the solvent to AT = 0.01K 
leads to a coarsening of the structure (Fig.[2ji). This 
coarsening continues as a function of time, when the 
sample is kept at the same temperature. This is demon- 
strated by Fig. [2^ which shows the same situation after 
a waiting time of a 20h. The process of the phase sep- 
aration taking place in Figs. [2^,-e can be followed more 
quantitatively by determining the number and type of 
the nearest neighbors for each particle, which we define 
as the number of particles with centers a distance < 1.1(7; 
away from a given particle of the same type (A,B). In 
Fig. [3] we display the fraction of each type of particles 
with a coordination number c > 2. We see that the 
fraction is low when the particles are homogeneously dis- 
tributed (Figsj2^-b), but grows rapidly to a value close 
to 1 when the system is phase separated due to the fact 
that in this situation it is more likely for particles to be 
surrounded by other particles of the same species (Figs. 
[2]i-e) . The coordination number c — 2 was the most fre- 
quently found configuration in the snapshot Figj2];, for 
both kinds of particles, with N c = 2 /N tot ranging from 0.3 
to 0.4. This coordination number mostly corresponds to 
string-like colloidal structures which we believe are the 
initial stages of phase separation, forming the beginnings 
of a bicontinuous network. 

For xa — 0.28 the corresponding experiments do not 
show the formation of a bicontinuous network and the 
CL-CL demixing instead proceeds via the growth of small 
clusters of the minority phase (A particles) which are sur- 
rounded by a liquid composed of the majority (B-type) 
particles. Similar to the case described above, the ar- 
rangement displayed in Fig.[2j (formed after « lh) hardly 
changes when the sample is kept at the same tempera- 
ture for a further 20h (Fig.[2]j) as the coarsening is rather 
slow. Such a strong dependence of the morphology of the 
demixing process on the particle composition has been 
already observed by other authors [22]. Here, however, 
the CL-CL demixing transition and associated critical 
point is induced by another critical point (that of the 
solvent), which demonstrates a coupling of critical phe- 
nomena over many decades in length scales. Additional 




FIG. 3: Fraction of particles N c> 2/N to t with coordination 
number c > 2, where Ntot is the total number of particles, 
calculated separately for each kind of particles (AA and BB), 
for the snapshots in Figs. [2k-e. 



evidence for the nature of the CL-CL demixing process 
may be obtained from Fig. [4j where we show snapshots 
for various different compositions at almost identical val- 
ues of AT. We observe a change from the 'island' mor- 
phology at lower values of xa to a bicontinuous structure 
when xa ~ 0.5. 

Having established from the experimental images (e.g. 
by comparing Figs. [2^, and[2ji) that approaching the sol- 
vent critical temperature leads to colloidal phase separa- 
tion, we now develop a model for the system free energy 
in order to further understand the phase behavior of the 
system and to construct a phase diagram. For simplicity, 
we consider a mono-disperse, two-dimensional systems 
with the following particle interactions: 



Uij{r) 



r < Veff 

r > a eff , 



(2) 



where the indices i,j = A,B and &jj(r) is the corre- 
sponding critical Casimir potential JI]) acting between 
the colloidal particles (we use the universal scaling func- 
tions &ij(z/£) from Ref. [6j) and the electrostatic repul- 
sion is modelled via a hard disk repulsion with effective 
diameter <J e ff = 3.1/im. For the other parameters we use 
the corresponding experimental values (see above). 

We construct an approximation for the Helmholtz free 
energy of the system by splitting the pair potential into 
the hard-disk part Uhd(f) plus the tail term: 



Wij{r) 




r < Veff 

r > a eff . 



(3) 



Thus, the full pair potential between the particles is 
u ij( r ) — u hd{r) + Wij(r). We assume the following 
mean-field (Van der Waals like) approximation for the 
Helmholtz free energy per unit area: 

!{PA,Pb) = fhd(pA, PB) + ^ PtPjWiJi ( 4 ) 
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FIG. 4: Snapshots of binary colloidal systems of three different compositions, xa = 0.28, xa = 0.32 and xa = 0.54, obtained 
after lh, at temperatures close to T c , with AT =0.01K. 
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FIG. 5: Phase diagram for the system plotted as a function 
of AT versus concentration xa = Pa/(pa + Pb), for colloids 
with a total density Pao\ + Pbo\ = 0.7. The index A and 
B correspond to the (+) and (— ) particles in the snapshots 
displayed in Figs. [2] and [4] We see that the theory predicts 
that for AT < 0.067K the Casimir forces between the colloids 
are strong enough to induce phase separation. 



where ibij — J dr u>ij (r) and fhd is the Hclmholtz free 
energy per unit area of a uniform fluid of hard-disks with 
bulk densities pa and ps ■ We have used the scaled parti- 
cle approximation for fhd - see Ref . [23] . fhd includes the 
(exact) ideal gas contribution fu = ksTJ2i pi(\n[A 2 pi] — 
1), where fee is Boltzmann's constant and A is the ther- 
mal de Broglie wavelength. 

As the temperature of the solvent approaches T — > 
T c . the strength of the Casimir potential increases and 
eventually leads to a CL-CL demixing transition. Within 
the simple model described above, the phase behavior 
depends on the values of the integrals Wij , and therefore 
on AT. In Fig. [5] we display the resulting phase diagram 
for the case when the total density pa&a + Pb&b = 0.7, 
which corresponds to the experimental results in Figs. 
[2] and [4] Because the strength of the critical Casimir 
forces are identical among like species of particles, this 
results in a symmetric demixing line (binodal) centered 
on the mixing ratio xa = 0.5. We see from Fig. [5] that 
as AT is decreased, if the mixture has the composition 
xa = 0.5, then one hits the critical point for the CL- 
CL demixing transition at AT ss 0.067A". We also see 



that the onset of demixing for a range of values of xa 
centered on xa — 0.5 occurs at a similar value of AT. 
This is in fairly good agreement with the experiments for 
xa = 0.54: looking at Fig. [2] one would estimate that c), 
which is for AT = 0.08K, is near the onset of demixing. 

The position of the critical point in Fig. [5] is fairly 
sensitive to the size of the electrostatic contribution to 
a e f f . Bearing in mind also some of the simplifications in 
our model, e.g. that we assumed the same size of both 
particle types, the difference between the experimentally 
determined CL-CL demixing point from the theoretically 
calculated value is surprisingly small. Therefore, it is 
clear that the results in Fig. [2J;, whilst not being exactly 
at the critical value of AT, are certainly for a state point 
nearby. 

One issue that we should comment on further is that, 
as mentioned above, the demixed particle configurations 
of our two-dimensional fluid change very slowly over time 
- compare, for example, the snapshots in Fig. [2ji and e, 
which are taken after lh and 20h, respectively. This is 
due to that fact that the characteristic length scale L(t) 
of the domains is a power-law function of time L(t) ~ t n , 
where we expect the exponent n « 1/3 (2U |5S]. This 
means that to observe further coarsening in the size of 
the domains, one would have to wait over much longer 
time periods. Furthermore, the dynamics is also slow 
due to the critical slowing down that one observes in 
the vicinity of a critical point. We have used dynamical 
density functional theory (DDFT) [57] , based on the 
free energy functional which for a uniform fluid gives the 
free energy in Eq. Q, in order to study the demixing 
dynamics as the binary colloidal fluid is quenched into 
the spinodal region of Fig. [5] We see qualitative agree- 
ment between the DDFT results (not displayed) and the 
snapshots from the experiments, displayed in Figs. [2] and 
[J] However, matching the DDFT to the experimental 
snapshots is not straightforward because an input to the 
DDFT is the long time diffusion coefficients Da and Db 
of the particles over the surface. One can estimate L(t) 
from a time sequence of experimental images, allowing us 
to estimate Da ~ Db = D and then one finds that the 
DDFT density profiles have roughly the same L(t) and 
appear qualitatively similar to the experimental snap- 
shots. 

In this Letter we have presented experimental and the- 
oretical results demonstrating that one can use critical 
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Casimir forces to precisely control the phase behavior, 
structure and dynamics of binary mixtures of colloidal 
particles. In particular, we have observed that by taking 
the solvent temperature close to T c , one can generate a 
further critical point in the binary colloidal mixture. This 
second critical point and the associated CL-CL dcmixing 
transition is generated entirely by the solvent critical fluc- 
tuations and resulting Casimir forces due to the tempera- 
ture proximity to T c . This striking observation also leads 
one to speculate that it would be interesting to make a 
system to which are added two additional species of much 
larger colloids. The Casimir force between these larger 
colloids, due to the critical fluctuations of the smaller 
colloids, could then induce a third critical point in the 
system, although, of course, depletion forces would also 
be important in such a system. 
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